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Abstract 

Background: The diffusion-limited reaction rate of a uniform spherical reactant is generalized to anisotropic reactivity. 
Previous work has shown that the protein model of a uniform sphere is unsatisfactory in many cases. Competition of 
ligands binding to two active sites, on a spherical enzyme or cell is studied analytically. 

Results: The reaction rate constant is given for two sites at opposite ends of the species of interest. This is compared 
with twice the reaction rate for a single site. It is found that the competition between sites lowers the reaction rate 
over what is expected for two sites individually. Competition between sites does not show up, until the site half angle 
is greater than 30 degrees. 

Conclusions: Competition between sites is negligible until the site size becomes large. The competitive effect grows 
as theta becomes large. The maximum effect is given for theta = pi/2. 



Background 

The purpose of this article is to generalize the Smolu- 
chowski [1] calculation of the steady state bimolecular 
rate constant, 4ttDR, to anisotropic reactivity. This was 
necessary because previous work [2,3], has shown that 
diffusion-controlled protein-ligand binding is not mod- 
eled successfully by a uniformly reactive protein for 
some cases. As a consequence of this, more realistic 
models of proteins were created. Previously this was 
done by having one reactive site on a spherical molecule 
or cell [4-6]. Berg and Purcell considered many sites on 
a sphere [7]. This paper considers two axially symmetric 
sites on opposite ends of a sphere reacting with small 
molecules. See Figure 1. This might be used to model a 
protein-ligand reaction: such as the diffusion-controlled 
reaction of acetylcholine with the active sites of tetra- 
meric mouse acetylcholinesterase [8]. 

Methods 

We want to solve 

3 2 c(r, Q)/dr 2 + (2/r)3c(r, 8)/3r 

+ l/r 2 (l/sin Q)d/dQ(sin Qdc (r, 8)/98) 

= 0 (1) 
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where c(r, 6) gives the ligand concentration. Figure 1 
shows that there is symmetry between the two active 
sites. As a consequence of this, the net flux normal to 
the equatorial plane vanishes for all r. This reduces the 
two site problem to a one site problem, where now 6 
varies from 0 to tt/2. Mathematically this says 

ac(r,e)/ae = o e = n/2 (2) 

We also have 

c 0 = lim r ^ooc(r, 6) (3) 

where c o is the bulk concentration of ligand. The reactive 
boundary conditions are given by 

c(R, 9) = 0 0<9<6, (4) 

and 

dc(r^)/dr\ r=R = 0 6o<e<7r/2 

(5) 

As done previously [6,9], equation (4) is replaced by 
the constant flux boundary condition 
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Figure 1 A non-reactive sphere with two reactive spherical caps 
of half angle 0° at each end. 



dc(r,Q)/dr\ r=R = Q 



o<e<e. 



(6) 



where Q is evaluated, by requiring the average concen- 
tration of ligand to vanish over the active site, (see ref. 6, 
equation 8, k — > <*>) 



Jo 



c(R, fysinddQ = 0 



(7) 



The diffusion limited rate constant for the two site 
problem is 



KDC 



2jtR 2 D r 
c° Jo 



-I 



=R sinB^e (8) 



Where D is the diffusion coefficient of the ligands and R 
is the radius of the sphere. The accuracy of the constant 
flux boundary condition may be seen, for small 6 o (e.g. 
small binding sites), by considering a reactive disk in an insu- 
lating plane. For this problem, the exact solution is known 
[6]. It is k DC = 4 Da, where a = the disk radius. The constant 
flux method yields [6] k DC = 3.7 Da. Thus the constant flux 
boundary condition is accurate for small reactive sites. 

The solution of equation (1) is given by 



c(r, 6) = a + Y7 m a ™fm( r ) p m{ cos9) 



(9) 



where the P m (x) are Legendre polynomials of order m. 
Application of the boundary condition given by equation 
(2) yields 



c(r, 6) = a + ^2 m a 2m f 2 m( r ) p 2m(cosQ) 



(10) 



Equation (3) is satisfied with a = c o and with the radial 
functions 



/ 2 » = l/r 



,2m+l 



(ii) 



The coefficients a 2m in equation (10) are found using 
equations (5) and (6). Details, parallel a previous deriv- 
ation [6]. 



Diffusion Limited Rate Constant vs. Active Site Half Angle 
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Figure 2 Plot showing the variation of the diffusion-limited rate constant with half-angle theta for 3 different cases. 
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Results and discussion 

In this section we present the results and discuss how 
the method represents a step forward in the field of 
diffusion-controlled kinetics and biophysics in general 

The results are given by the solution of the model in 
the preceding section, and by comparing it to twice the 
rate limited constant for the one site problem [6]. The 
difference between the two problems, gives a measure 
of competition between the two sites for reaction with 
ligands. 

The diffusion limited rate constant for the two site 
problem is 



k2 = k D c/4nDR 



(12) 



where 4ttDR is the rate constant for a uniformly reacting 
sphere. The reciprocal of our rate constant is 



l/k2 = 1 



4(1- cos e„) 2 

(P 2m -i ( cos 6.) -P 2m +i ( cos 6„)) 2 / (m + 0 (im + ^ 



(13) 

The reciprocal of the rate constant for one site is given 

1 v" 

2fl-cos0„f ^ m=1 



by 

1/kl = 1 + 



[P m .! ( cos 6.) -P m+1 ( cos e.)) 2 / (m + ^ (m + 1) 

(14) 

Figure 2 shows plots of I<2 (8 o ) (the rate constant for 2 
sites), twokl (6 o ) (twice the rate constant for 1 site) and 
kl (6 o ) (the rate constant for one site [6]) versus 6°. For 
small 6 o , the two site rate constant and 2 x the single site 
rate constant are in agreement with each other, as would 
be expected(both behaving as two active sites on a large 
sphere). As 6 o grows above 30 degrees, the curves grow 
apart. The two site curve, being less than 2 x the one site 
curve. 

The difference between the curves is a measure of the 
competition between the two sites. The competition effect 
does not show up till around 30 degrees. For 0 o = f, 
k2 = 1, which is the exact result. Thus, the constant flux 
boundary condition, is good for large 0 o . 

This paper represents progress in the field, by present- 
ing a new model for the interpretation of experimental 
data. This is for macromolecule-ligand binding reactions 
that fall in the diffusion-controlled regime. Previously, 
only one site models were available for modeling pro- 
teins. Proteins with multiple binding sites can now be 
studied. 



Conclusions 

The analytical expression for the diffusion-limited rate 
constant to two active sites on a sphere has been given. 
The result was used to study the competitive effects 
between the two sites. The effect doesn't show up until 
the site half-angles is greater than 30 degrees. The 
competitive effect grows until its maximum value is 
reached at 6 o = n/2. 
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